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We test the cooling algorithm with gluonic and staggered hadronic spectroscopy on SU (3) gauge field configu- 
rations generated with two flavors of staggered dynamical fermions. We find cooling is not reliable as the basis for 
improved hadronic operators. We also find that performing cooling sweeps to reveal more clearly the topological 
properties of the gauge fields eliminates the spin structure of the hadron spectrum. 



1. Introduction 

Motivated by Shuryak et. al. [1] who have ar- 
gued that instantons are the driving mechanism 
for particle mass generation in QCD, Chu, et. al. 
[2] used the cooling algorithm (p = co Metropo- 
lis) on quenched SU (3) configurations and mea- 
sured hadronic spatial correlation functions and 
various gluonic quantities. They performed 0, 25 
and 50 cooling sweeps. They claim to find rea- 
sonable agreement with the work of Shuryak et. 
al. for cooled and uncooled results, and argue 
in support of the dominant role of instantons. 
However, Trottier and Woloshyn [3] calculated 
hadronic spectroscopy on cooled SU (2) config- 
urations (which converge to a global minimum 
far more quickly) and large m^a. They find the 
mesonic spectrum becomes degenerate after a few 
cooling sweeps. 

Why use cooling? Operators that measure 
topological properties are plagued by ultravio- 
let (short distance) effects. Cooling eliminates 
these effects and exposes the metastable states 
of instantons, but the method has bias! Infinite 
cooling recovers a free field configuration. Cool- 
ing, or rather smearing, in dimensions orthogonal 
to time is valid for glueball states, but not for 
hadronic states. Teper [4] has argued that cool- 
ing is not uniform and short distance modes re- 
lax more quickly than long distance modes (which 
suffer from critical slowing down) to a global min- 
imum. Therefore, we might expect to see short 
distance measurements converging to free field 
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values faster than long distance measurements. 

In this paper we will view the cooled gauge 
fields as composite fields (defined through the 
cooling algorithm!). We will test the hypothesis 
that hadronic and gluonic operators defined on 
these composite fields are "improved" operators: 
does the statistical improvement from (a few) 
cooling sweeps outway the inherent bias while at 
the same time exposing the long distance physics, 
namely instantons? Detailed results will be pre- 
sented elsewhere [5]. 

For our tests, we used 200 of the /3 = 5.6, 
niqa = 0.01 & 0.025 2-fiavor staggered fermion 
gauge field configurations at SCRI generated by 
the HEMCGC collaboration. We performed ex- 
treme amounts of cooling sweeps, namely 0, 5, 
25, 50, 100 and 500, to study its effects. For spec- 
troscopy, we used (four fiavor) staggered fermions 
with valence quark mass equal to the sea quark 
mass. We have also measured the eigenval- 
ues of the Dirac operator (on selected configu- 
rations), {tptp), glueballs and Polyakov loops, and 
the heavy quark potential. All calculations were 
performed on the CM-2 at SCRI. 

2. Simulations 

In figures 1 and 2 we show the heavy quark 
potential V(R)a calculated at time slices T = 4, 
5, 8, 10, 12 for 5 and 25 cooling sweeps. We de- 
termine the lattice spacing by the force method 
of Sommer and compare to Ref. [6]. We find the 
string tension and lattice spacing are unchanged 
within errors between and 5 cooling sweeps. 
However, after 25 cooling sweeps we see a clear 
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/S = 5.8. am = 0.01. T=4, 5, 8,10.12 cool=5 
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Figure 1. Heavy quark potential V(R)a after 5 
cooling sweeps. The potential is shown at T = 4 
(o), T = 5 (x), T = 8 (+), T = 10 (O), T = 12 
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Figure 2. Heavy quark potential V(R)a after 25 
cooling sweeps. Symbols are the same as in Fig. 1. 
Note the separation of V(R)a at smaller T. 



0.15 



e = 5,6, m„a = 0,01, cool=500 



-Q(t) 
-s(t) 

- <7T(t)n(0)> 
-<7T(t+16)7T(16)> . 



0.10 




0.05 — / 



0.00 



10 



20 



30 



Figure 3. Time profile of Q(t), S(t), and the pion 
propagator (7r(t)7r(0)) and (7r(t + 16)7r(16)). 
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Hadron masses at m„a 



0.01. The confidence 



levels are not shown but are greater than 68%. 
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separation of the potential from T = 4 and 5 but 
unchanged at T >= 8 except for an overall scale 
change - consistent with cooling altering short 
distance modes faster than long distance modes 
(non-uniform minimization). 

We also observed this phenomena in our mea- 
surements of the 0"*""*", 1"' , 2+"'" glueballs masses 
using different levels of 3d smearing in addition to 
cooling. These operators are variational bound 
states, and if we had not smeared in the time 
direction we would expect to see the glueball ef- 
fective mass bound decrease in time. However, at 
fixed cooling sweeps, we find the glueball effective 
mass increase in time, and at fixed time the effec- 
tive mass decreases under cooling consistent again 
with non-uniform minimization. Comparing to 
Ref. [8], att = 2, meff(0++) = 0.82(5), 0.46(2), 
0.29(2), 0.24(2), 0.20(1), 0.118(8). Only at 500 
cooling sweeps can we reliably extract m(0"'""'") = 
0.23(4)[0.40/12]. _ 

We calculated {tptp) using one Gaussian field. 
We found (V^V) = 0.1122(5), 0.0202(4), 0.0283(6), 
0.0278(7), _0. 0265(9), 0.0249(8) for m^a = 
0.01, and (V't/') =0.2140(3), 0.0589(2), 0.0541(2), 
0.0523(2), 0.0508(2), 0.0491(2) fo£ m^a = 0.025 
as a function of cooling. We see {tptp) is decreas- 
ing but not as fast as in SU (2) [3] - we see no 
evidence that the instanton induced zero modes 
for staggered fermions are not becoming delocal- 
ized as required for chiral symmetry breaking [7] . 

In Fig. 3 we show a time profile of the topolog- 
ical charge density and energy density summed 
over 3-space, along with the wall source point 
sink zero momentum pion propagator with wall 
source at t = and 16 on a single configuration 
at 500 cooling sweeps. We see the charge and 
energy density (normalized in units of a single in- 
stanton) significantly overlap and correspond to 
large bumps in the pion propagator. Therefore, to 
enhance statistics and restoration of translation 
invariance we measured our spectroscopy with 
16 time-slice sources. We show the results for 
m^a = 0.01 and 0.025 in Tables 1 and 2. We 
see the mesons for both m^a are becoming de- 
generate indicating that cooling has eliminated 
spin structure. In addition, fiavor symmetry is 
restored after 5 cooling sweeps. The nucleon is 
consistently about 3/2 times larger than the rho. 



3. Conclusions 

We find cooling is not reliable as an improved 
hadronic operator. A few cooling sweeps leaves 
long distance gluonic physics intact (as is well 
known) but has dramatic effects on the staggered 
hadron spectrum. As in the results of ref. [3], 
we find that performing more cooling sweeps to 
reveal more clearly the topological properties of 
the gauge fields eliminates the spin structure of 
the hadron spectrum. As a function of cooling we 
find the pion mass makes a precipitous drop (in- 
dicating heavy mixing with short range fiuctua- 
tions) then slowly rises. All other particle masses 
decrease uniformly and the overall spectrum be- 
comes degenerate. We did not perform enough 
cooling sweeps for the spectrum to become free 
field like. 

It is not clear what has been revealed about the 
instanton content of the QCD vacuum. Instan- 
tons definitely affect quark propagators as can be 
seen from the eigenvectors of the dirac operator 
or the profiles of the pion propagator (Fig. 3). 
This is not too surprising since very little energy 
content is left. However, cooling is definitely elim- 
inating staggered fermion spin structure which is 
undoubtedly related to the large energy remain- 
ing after a few cooling sweeps. Without an a pri- 
ori cutoff, the utility of the cooling algorithm is 
questionable. 
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